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Abstract 

Using the self-consistent model, we present nonlocal spin-transfer effects caused by the feedback 
between inhomogeneous magnetization and spin-transfer torque on the current-induced magneti- 
zation dynamics in nanomagnets. The nonlocal effects can substantially improve the coherence 

i_> 

time of precession in nanomagnets and thus reduce the linewidth of power spectrum. This narrow 
linewidth results from the nonlinear damping of spin-waves due to the nonlocal spin torque which 
is inherent and thus should be considered in future experiments. 
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A spin-polarized current exerts a torque to a ferromagnet (FM) by transferring the spin- 
angular momentum, i.e. spin-transfer torque [l| (STT). STT generates a new class of mag- 
netization (M) dynamics in spin- valve structures 0, y] , potentially useful for applications 
in magnetic nonvolatile memories and microwave oscillators. 

In layered structures where the current flows perpendicular to the plane, the direction and 
magnitude of STT at a point r is decided by the spin accumulation /x s and associated spin 
current 3 S at the same point r. Initial theories [1] assumed that the dependence of /x s on M is 
local and thus essentially fixed by the local M at the same point r. However, its dependence 
on M is inherently nonlocal because of the 3-dimensional (3D) spin diffusion 4, ^, 6]. In 
other words, when the conduction electron arrives at a point r on the FM| normal metal (NM) 
interface, the reflected (transmitted) electron takes the spin direction anti-parallel (parallel) 
to the local M at the point r, diffuses along the interface, and then transfers its spin-angular 
momentum to another local M at a far away point from the r. That is, /j, s at a point r is 
affected by all local M's at other points. The local assumption becomes really invalid when 
M is inhomogeneous. Note that micromagnetic [l,!^ and time-resolved imaging studies 9] 
have revealed excitations of incoherent spin-waves and thus inhomogeneous M due to STT. 
In this situation, the effect of ii s on M (=STT) and the nonlocal effect of M on the ii s should 
be treated on an equal footing. The conventional treatments, which ignore the latter part, 
actually deal with only half of the relevant parts. Therefore, the self-consistent feedback 
between inhomogenous M and STT through the nonlocal effect should be considered. 

The STT caused by the nonlocal effect can be named nonlocal spin-transfer torque 
(NLST) since it allows a single FM with inhomogeneous M to exert spin-transfer effects 
on itself. Despite efforts to investigate NLST, the understanding of the M dynamics af- 
fected by NLST remains elusive especially for the spin-valve structure which is important 
from the viewpoints of fundamental physics and applications. Previous theoretical studies 
on NLST 4j, |5|, |6j have addressed the phenomenon in the perturbative regime of small spin- 
wave amplitudes and thus could not investigate the dynamic modes for the current exceeding 
the threshold for the onset of magnetic excitation. Previous numerical studies 10] lacked 
the exact calculation of 3D dynamic motion of n s and focused only on the single FM. 

In this Letter, we have directly calculated 3D dynamic motion of fi s self-consistently 
coupled with the M dynamics, which allows us to apply the model to both single FM and 
spin-valve structures. This self-consistent treatment is essential to correctly describe unique 
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spin-wave modes caused by NLST and explains two important experimental results: spin- 



wave excitations in a single FM 11] and narrower linewidths in spin- valves than are expected 

□ 

within the assumption of homogeneous M [12j • 

The equations of motion of M (Eq. (1)) and fj, s (Eq. (2)) are self-consistently 

solved for FM and NM. 

d t m = — 7_f(iti x H e jj) + am x d t m (1) 
+ j F /(M s t F ) J s | ~t F /2 — Js\+t F /2\ , 
dtfJls + V • J s = -7at(Ms X H exf ) - fls/Taf. (2) 

Here m is the unit vector of M, 7f(tv) is the gyromagnetic ratio of FM (NM), H e // is 
the effective field including the magnetostatic, exchange, external (H ext ), current- induced 
Oersted, and thermal fluctuation fields, a is the intrinsic damping constant, M s is the 
saturation magnetization, tp is the thickness of FM, J s = —DVn s is the spin current, D 
is the diffusion coefficient, r s / = l 2 s f/D is the spin-flip scattering time, and l s f is the spin- 
diffusion length. The change of charge and spin current ,L and J s at the interface of FM|NM 
are related to the potential drop over the interface as [l3j 

J e = (G T + G^AfiJe + (G T - GJm • (A/*./e) (3) 
J s = (h/2e 2 )[Re(G n )m x (m x 2Afi s ± ^m), 

- ((G T + G x )m • A Ms - (G T - G;)Aju e ) m], (4) 

where /i e is the electric potential, A/i = /i(±^/2+0)— /x(±t j p/2— 0) is the potential drop over 
the interface, G s (s=| or {) is the spin-dependent conductivity, = (Gf — G^)/(Gf + Gj_) 



is the bulk (inter 

n 

is disregarded [14J |. 



ace) spin asymmetry, Gfj_ is the mixing conductivity. A small Jm(G^) 
At the interface of FMINM, J P and J, • m are continuous under the 



condition of /j, s x m = in FM. /j, s and m are related through the Eqs. (2)-(4), and the 
spin-version of the Ohm's law with the boundary conditions of \i e = — eV(0) and /x s = 0(0) 
at the far-right (-left) end of the NM electrodes. 

To validate the self-consistent model, we first carried out simulations for the single FM, 
Cui(10)|Co(tc , o)|Cu 2 (52 — t) (all in nm) where tc Q varies from 2 to 8nm, and compared 



modeling results to the experimental ones in the Ref . 11[ . Since this structure has no second 
FM, the conventional LLG-Slonczewski equation is not applicable. Asymmetric Cu leads 
provide asymmetric /x s at each side of the Co layer (Fig. 1(a)). Jl at interfaces (= ^ U1 ^ c ° + 



^Co\Cu 2 ^ j g ne g a ^j ve w h en the electron flows from the thick to thin Cu layers, corresponding 
to a negative current. This negative JI provides negative NLST. Fig. 1(b) shows the time 
evolution of averaged out-of-plane component of M ((M z )) at various negative currents when 
the out-of-plane field H is 2.5T. M initially saturates along the out-of-plane direction, but 
cannot keep the saturation state at negatively large currents even when H is larger than the 
out-of-plane demagnetization field (w 1.6T) (Fig. 1(c)). When the current is turned on, 
a tiny in-plane component of M is developed especially at the long edges where the Oersted 
field is the largest. Interplay between this inhomogeneous M and negative NLST excites 



spin-waves, resulting in t 
As in the experiment 



re rapid decrease of (M z ) within a few nanoseconds. 



Ill ] . we observed current-induced excitations only at negative cur- 



rents. The normalized modulus of the magnetic moment |M| is much smaller than M s at 
those bias conditions (Fig. 1(d)), indicating excitation of large amplitude incoherent spin- 
waves. Inset of Fig. 1(d) shows a snapshot of domain pattern at H = 2.5T and I = —11mA. 
Local M's at edges are mostly in the plane whereas those near the center of cell are in vortex- 
like states caused by negative NLST which prefers non-collinear configuration of local M's. 
This inhomogeneous M results in the reduction of the average spin accumulation in the NM 
and thus the reduction of resistance of the stack (not shown). 

When H > H^, the critical current Iq for excitations linearly depends on H (Fig. 1(c) 
and (d)). As shown in Fig. 1(e), numerical results of the slope (= dlc/dH) are in better 
agreement with the experimental ones than the theoretical ones (for the theoretical Iq, see 
Eq. (10) in the Ref. Q]). In the experiment fll| . the intercept of extrapolated boundary 
at / = is nearly zero for the sample with tco — 8nm, whereas the theoretical intercept is 
about 0.8T(m Hd/2) for all thicknesses (inset of Fig. 1(e)). For tc = 8nm, the numerical 
intercept is 0.23T and again in better agreement with the experimental one. We attribute 
these better agreements to the fact that the self-consistent model more realistically takes 
into account the influence of the shape and finite size of nano-pillar on the spin-wave mode. 

Fig. 1(f) and insets show eigenmode analysis for the M dynamics at / = —11mA and H = 
2.5T, corresponding to a periodic oscillation of (M z ). As shown in Fig. 1(b), however, the 
(M z ) oscillation is in general nonperiodic for most negative currents due to highly nonlinear 
coupling among local M's through the NLST. A rule of the bias condition for a periodic 
oscillation may exist but we could not find it because of a fixed step size of / and H in our 
simulations. At this bias condition, the power spectrum shows two peaks at /l(= 75.3GHz) 
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and 2/x, where fi = JcoH '/2ir. The eigenmode images (insets) show that the precession 
region with a higher power is localized at edges. Note that these eigenmodes are unique 
features of the NLST and not expected in the field-driven excitation [161 ] . 

This result demonstrates one crucial implication of the NLST, namely, destabilization 
effect of negative NLST on local M's. Via the spin diffusion, the electrons backscattered from 
the FM destabilize local M's whereas the electrons transmitted through the FM stabilize. 
The two effects always exist simultaneously but one dominates the other because the NM 
electrodes and thus fj, s are not symmetric. Since the sign of /j, s is reversed by changing the 
current polarity, the stabilizing effect is expected for a positive current, i.e. positive NLST. 
In the single FM excitation, we observed almost macrospin behaviors for positive NLST. In 
a spin-valve, however, different types of spin-wave modes are expected for positive NLST 
because the local STT is nonzero and thus generates incoherent spin-waves |8j. 

In the second study, we applied the self-consistent model to a spin-valve structure, 
Cu(80)|Py(20)|Cu(6)|Py(2)|Cu(2) (Py=Permalloy) experimentally studied by Sankey et 
al. 12j. They reported a surprising result that the current-induced dynamic modes can 
generate narrower linewidths at low temperatures than those expected within the macrospin 
assumption. To investigate the origin of this experimental finding, we performed simulations 
with three different approaches: i) macrospin model (MACRO), ii) conventional micromag- 
netic model without considering NLST (CONV), and iii) self-consistent model (SELF). Fig. 
2(a), (b) and (c) show contours of spectral density of (M x ) as a function of / at AK when the 
effective field of 500Oe is applied along the in-plane easy axis (/ fx). The positive current 
corresponds to the electron-flow from Cu(2) to Cu(6), and thus positive NLST. MACRO 
shows the well-known red- and blue-shift depending on I (Fig. 2(a)). CONV shows only 
red-shift up to a critical current (Iq ONV « 2mA, Fig. 2(b)). When / > I^ ONV , M dynamics 
in CONV becomes complicated due to excitations of incoherent spin-waves. As indicated 
by an arrow, we observed secondary peaks with about half the frequency of main peaks, 
corresponding to the precession of end domains 8]. In SELF, we observed similar secondary 
peaks indicating non-single domain state, but much clearer peak structures than CONV up 
to about 2AmA which is larger than I^ ONV (Fig. 2(c)). It indicates that the positive NLST 
provides a more periodic oscillation than that obtained in CONV. 

Fig. 2(d) shows power spectra obtained in the three models. At a low temperature 
(T), SELF shows the narrowest linewidth whereas CONV produces the broadest one due 
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to excitations of incoherent spin-waves. We calculated the T dependence of linewidth from 
lorentzian fits (Fig. 2(e)). At low temperatures, SELF provides narrower linewidths than 



MACRO, consistent with the experimental observation [12|. Therefore, the positive NLST 
indeed results in a substantial improvement of the coherence time of precession although M 
is not in the single domain state. It indicates that it is possible to reduce the linewidth by 
properly controlling the NLST. In MACRO, the linewidth monotonously increases with T. 
On the other hand, in SELF, the linewidth linearly depends on T for T < 50K and more 
rapidly increases at higher temperatures. 

The narrower linewidths in SELF are caused by two nonlinear effects of the positive 
NLST: an increase of the effective exchange stiffness in short range and an increase of the 
damping of incoherent spin- waves in long range. As a result, the positive NLST provides 
an additional nonlinear spin- wave damping. For a spin-torque nano-oscillator, the linewidth 
Au in the low-temperature limit is given by 



= r + (p ) (k B r/E ) [i + (ivyr e// ) 2 ] (5) 

where N = du(P)/dP is the nonlinear frequency shift coefficient obtained from u>(P) = 
co>o + NP, co>o is the ferromagnetic resonance frequency at I = 0, P is the normalized power, 
T e ff — a (I + QI C ) is the effective nonlinear damping, Q is a phenomenological coefficient 
characterizing the nonlinear positive damping, and I r is the critical current for the magnetic 



excitation (for details of other parameters, see Ref. 18] ). 

Eq. ([5]) predicts two important consequences of the nonlinearity. First, the linewidth of 
an auto-oscillator with a nonlinear frequency shift (i.e. N ^ 0) increases by the factor (1 + 
(N/T e ff) 2 ) from that of a linear oscillator (i.e. N = 0). Second, the linewidth of a nonlinear 
oscillator decreases with increasing the nonlinearity of damping Q. It is because the linewidth 
is determined by nonlinear properties of the system where the normal linear damping is 
compensated by local STT. In this case, an increase of the nonlinearity of damping can lead 
to a decrease of the linewidth, known as the noise suppression due to nonlinear feedback 

Inset of Fig. 2(e) shows that N is nonzero and almost identical for the two models. 
Thus, the linewidth is wider than that expected in a linear oscillator. Using Eq. ([5]), we 
fit the values of Q from the calculated linewidths at T = 107^ and obtained Q = 0.13 in 
MACRO and Q = 1.96 in SELF. The fit value Q in SELF is consistent with the assumed 



19] 



values (Q — 1 ~ 3) [18| to explain experimental observations. Note that the nonlinear 
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theory referred here does not take into account the spin transport and the Q value has been 
used as a fitting parameter without justification of its origin. In contrast, our self-consistent 
treatment shows that the large Q is mainly caused by NLST. Thus, we conclude that the 
nonlinear spin-wave damping due to NLST is responsible for narrower linewidths in SELF 
at low temperatures. When T is too high, the thermal random force overcomes the nonlocal 
effect due to positive NLST and thus the linewidth abruptly increases. For the opposite 
current polarity (i.e. negative NLST), we observed an increase of the linewidth (not shown). 

Finally, we note that the magnitude of NLST is easily controlled by modifying the asym- 
metry of layer structure like the conventional local STT. The effect of NLST on the current- 
induced M dynamics is determined by the ratio of NLST to local STT. Fig. 2(f) shows 
the ratio at the parallel magnetic configuration as a function of the thickness of Cu spacer 
(tcu) for the spin-valve structure studied here. The ratio is about 0.1 at tc u = Qnm which 



is the case of the Ref. 



121 ] . Note that the effect of NLST on the current-induced M dy- 



namics is considerable although the ratio is only 0.1. Furthermore, this ratio increases with 
increasing tc u as shown in Fig. 2(f). Therefore, NLST should be considered in designing 
and interpreting future experiments. 
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FIG. 1: (Color online) Current-induced excitation of single ferromagnet {tc = 2nm). (a) Spin 
accumulation pattern, (b) Time evolution of {M z ) at H = 2.5T and various negative currents, 
(c) and (d) show contours of (M z ) and \M\ as a function of H and /, respectively. Inset of (d) 
shows a domain pattern obtained at H = 2.5T and I = —11mA. (e) Slope of the critical boundary 
(= dlc/dH) as a function of tco- Inset of (e) shows the intercept of the extrapolated boundary as 
a function of tco- (f) Power spectrum at H = 2.5T and I = —11mA. Insets of (f) show eigenmode 
images for the two peak frequencies. Model parameters: Elliptical shaped pillar with 60 x 30nm 2 , 
M s = lA20emu/cm 3 , the exchange stiffness constant A ex = 2 x W~ G erg/cm, a = 0.01, the unit 
cell size= 3nm, and the discretization thickness of Cu layer varies depending on the total thickness 



and is not larger than 5nm. For Cu and Co, the spin transport parameters 



151 ] are bulk resistivity 



15 



nm 



2 s- 1 )=41 and 1.7. 



p(uQcm)=0.6 and 7.5, /3=0 and 0.46, l sf (nm)=450 and 59, and D(xl0 
For the interface Co|Cu, the parameters are interfacial resistance AR*(mflum 2 ) =0.51, 7=0.77, 
interfacial spin memory loss 5=0.25, and Re(G n )(x lO^Q" 1 cm =o.o. 
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FIG. 2: (Color online) Current-induced excitation of spin-valve. Contours of spectral density 
of (M x ) at T = 4K obtained in (a) MACRO, (b) CONV, and (c) SELF, (d) Comparison of 
power spectra obtained in the three models at T = WK. (e) Linewidth as a function of the 
temperature. Inset of (e) shows the frequency versus the power normalized by \M\. (f) The ratio 
of NLST to local STT as a function of the thickness of Cu spacer {pcu)- F° r (d) and (e), we chose 
different bias conditions for each model to yield a similar precession angle since the linewidth is 
proportional to the precession angle [3 ; ]. However, the main finding is not altered even when we 
choose the same bias condition for all models. Model parameters: Elliptical shaped nanopillar with 
120 x 60nm 2 , M s = 645emu/cm 3 [l7|], A e 



1.3 x 10 erg /cm, a = 0.025 [17j, and the unit cell 
size=5nm. For Py, the spin transport parameters are p(^Qcm) = 25.5, /3 = 0.7, l s f(nm) = 5.5, 
and D(x 10 15 nm 2 ,s _1 ) = 1.7. For the interface Py|Cu, the parameters are AR* (mf2/um 2 )=0.97 ', 
7=0.77, 5=0.25, and i?e(Gfj_)(x 10 10 il _1 cm _2 )=6.0. Parameters of Py were provided by Cornell 
group. Re(Gn) of Py|Cu was determined to mimic the critical current in the Ref. [12|. The pinned 
layer M is fixed along the in-plane easy axis (no stray field from it). 
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